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Abstract 

FLASH plans to use a "third harmonic" (3.9 GHz) superconduct- 
ing cavity to compensate nonlinear distortions of the longitudinal 
phase space due to the sinusoidal curvature of the the cavity volt- 
age of the TESLA 1.3 GHz cavities. Higher order modes (HOMs) 
in the 3.9 GHz have a significant impact on the dynamics of the 
electron bunches in a long bunch train. Kicks due to dipole modes 
can be enhanced along the bunch train depending on the frequency 
and Q-value of the modes. The enhancement factor for a con- 
stant beam offset with respect to the cavity has been calculated. 
A simple Monte Carlo model of these effects, allowing for scatter 
in HOM frequencies due to manufacturing variances, has also been 
implemented and results for both FLASH and for an XFEL-like 
configuration are presented. 
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1 Introduction 



FLASH plans to use a "third harmonic" (3.9 GHz) superconducting cavity to 
compensate nonhnear distortions of the longitudinal phase space due to the sinu- 
soidal curvature of the the cavity voltage of the TESLA 1.3 GHz cavities. Higher 
order modes (HOMs) in the 3.9 GHz have a significant impact on the dynamics of 
the electron bunches in a long bunch train. The analysis here seeks to determine 
what level of damping, if any, is required. 

In the case where the spacing of the bunches is near (but not exactly the 
same as) a multiple of the period of a dipole HOM, kicks from the HOM can 
be resonantly enhanced along the bunch train, depending on the frequency and 
Q-value of the modes. The enhancement factor for a constant beam offset with 
respect to the cavity can be expressed in a simple analytic form, which also 
provides the energy loss from resonance effects with a monopole mode. 

A simple Monte Carlo model of these effects, allowing for scatter in HOM fre- 
quencies due to manufacturing variances, has also been implemented and results 
for both a FLASH-like and for an XFEL-like configuration are presented. 

The beam parameters for our analyses are: 



Parameter 


FLASH 


XFEL-like injector 


Bunch spacing 


IfIS 


200ns 


Bunch charge 


InC 


InC 


Bunch length (la) 


Ips 


Ips 


Beam energy 


130 MeV 


500 MeV 


Bunch offset at entry 


1mm 


1mm 


Bunches per train 


800 


800 


Number of cavities 


4 


32 



Table 1: Beam parameters for typical FLASH and XFEL injector running. 

An important issue is the specification for how much kick can be tolerated 
before lasing stops. At FLASH, the spot size at the location where the "third 
harmonic" cavities will be installed is about 0.2mm and the invariant emittance's 
design value is 1 mm-mrad in both x and y; in real operation it is often twice that. 
The beam divergence is thus about 20/irad in the best case. We use ztlO/xrad 
as our target. We do not have a target for energy loss, but find small values for 
these effects in all cases. 

Although we are primarily concerned with FLASH, our methods are com- 
pletely general and we have done some investigation of the situation for XFEL. 
The optics for the XFEL are in a state of flux at this writing; our selection of 
parameters here is perforce somewhat arbitrary. The ±10/irad requirement is 
not far from other parameter sets that are under consideration at this writing. A 
more detailed study of the XFEL requirements is being undertaken by Yauhen 
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Kot and Thorsten Limberg. We do not here allow for betatron phase advance 
between the cavities, and this effect will be larger at the XFEL than at FLASH. 



2 Wakefield due to HOMs 

The purpose of this section is to define the parameters that are important for the 
long range wake field calculation. Our development follows reference [1] closely. 



2.1 Modes in a cavity 

2.1.1 The electric and magnetic fields 

Consider a monopole (m = 0) or dipole mode (m = 1) mode with the frequency 
/ = uj/{2tt) in a cavity with cylindrical symmetry. One obtains in complex 
notation for the electric and magnetic field: 

E{r,(j),z,t) = (^Er"^\r,z) cos(m{/)) Cr 
+ E^^\r, z) sin(m0) 
+ E^\r,z) cos(m0) ) exp(— zcjt) 



(1) 



B{r,(j),z,t) = \^Bi^"'\r,z) sin(m</)) e^, 

+ Bj^\r,z) cos(m0) 

+ Bi"^\r,z) sin(m0) ^ exp{—iut). 

2.1.2 The loss parameter and R/Q 

The interaction of the beam with a cavity mode is characterized by the loss 
parameter k^l^\r) or by the quantity R/Q |2]. These parameters can be deter- 
minated from the numerically calculated fields using the MAFIA post-processor 
[3l H] . The longitudinal voltage for a given mode at a fixed radius r is defined as 

Vi'^\r)= f"^ dzEi'^\r,z)exp{-tujz/c), (2) 

" JO 

while the total stored energy is given by: 
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(3) 
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From the voltage and stored energy the loss parameter and R/Q can be cal- 
culated: 



(4) 



For monopole modes the superscript (0) is usually omitted R/Q = R^'^^/Q. Al- 
though our definitions include a radial dependence in the loss parameter, R^^^ / Q 
is independent of the radius r since it can be shown (see [21 [5]) that V^"^\r) ~ 
and therefore k^"^\r) ~ r^™. 



2.1.3 The geometry parameter Gi and the Q-value 

The power Pgur dissipated into the cavity wall due to the surface resistivity Rsur 
can be calculated from the tangential magnetic field: 

Rsur 2 -^sur J dA . (5) 

For a superconducting cavity the surface resistance is the sum of the BCS (Bardeen, 
Cooper, Schrieffer) resistance Rbcsi which depends on the frequency and the 
temperature, and a residual resistivity Rq. The BCS resistance Rbcs scales with 
the square of the frequency / and exponentially with the temperature T: 

RBCs{f,T) oc ^ exp(-1.76r,/T). (6) 

The less- well understood residual resistance Rq adds directly to Rbcs but remains 
in the limit T — > 0. 

The total damping of a cavity mode is not only determined by the surface 
losses but also by coupling to external waveguides (HOM-dampers). Therefore 
one has to distinguish the Q-value Qo which is defined above and the external 
Q-value Qext which characterizes the coupling to external waveguides. Typically, 

Qo » Qext- 

The geometry parameter Gi [6] is defined as: 

Gl = Rsur Qo- (7) 

Gi is a purely geometric quantity that is independent of the cavity material; it 
depends only on the mode and the shape of the cavity that creates that mode. 
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Figure 1: A point charge qi traversing a cavity with an offset ri followed by a 
test charge q2 with offset r2 . 

2.2 Wakefields 

2.2.1 Wake potential 

Consider the situation shown in Figure [TJ A test charge q2 follows a point charge 
qi at a distant s. The distant s is positive in the direction opposite to the motion 
of the point charge qi. Both charges are relativistic (f ^ c). The Lorentz force 
on the test charge due to the fields generated by the point charge qi is 

F = ^= q2{E+cez x B). (8) 
The wake potential of the point charge gi is defined as: 

1 

W{x2,y2,xi,yi,s) = — / dz{E + c x B)t=(z+s)/c- (9) 

qi Jo 

The wake potential is the integrated Lorentz force on a test charge. Causality 
requires W (s) = for s < 0. 

The longitudinal and transverse components of the wake potential are con- 
nected by the Panofsky-Wenzel theorem [7] 

d 

±{x2,y2,xi,yi,s) = -V ±2W\\{x2,y2,xi,yi,s). (10) 

Integration of the transverse gradient (applied to the transverse coordinates of 
the test charge) of the longitudinal wake potential yields the transverse wake 
potential. 

2.2.2 Multipole expansion of the wake potential 

If the structure traversed by the bunch is cylindrically symmetric then a multipole 
expansion can be used to describe the wake potential. The location of the bunch 
train in the (r, 0) plane will break the symmetry and determine the azimuthal 
orientation of the m > modes. Consider again the situation shown in Figure [TJ 
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Assume that the point charge qi traverses the cavity at position (ri,(/?i), while 
the test charge follows at position (r2,(f2)- The longitudinal wake potential may 
be expanded in multipoles: 



Py||(ri,r2,</Pi,</^2,s) = J2 ^i"^2"^ ^ir(^) cosm(^2-V'i). (11) 

m=0 

The functions iy||'"'*(s) are the longitudinal m-pole wake potentials. There is no 
a-priori relation Between the wake potentials of different azimuthal order m. 

The transverse wake potential can be calculated using the Panofsky-Wenzel 
theorem, and the transverse m-pole wake potentials are defined as: 

W'r\s) = - r ds'wl.'^\s'), (12) 

for m > 0. There is no transverse monopole wake potential. The dipole wake 
potential does not depend on the position of the test charge q2- The kick on the 
test charge is linear in the offset of the point charge qi. 

2.2.3 Wakefields due to HOMs 

It is possible to write the m-pole wake potentials iy||™''(s) as a sum over all modes: 
Wf\s) = -Y. ^ni—g-j cos(t^„ s/c) exp(-l/r„ s/c) 

(13) 

—Q-j sin(cj„s/c) exp(-l/r„ s/c). 

where Un are the frequencies of the m-pole modes. A damping term has been 
included with the damping time r„ for mode n. As Qo >> Qext, the damping 
time of the voltage is very nearly 

r„ ^ ll^£^!k. (14) 



3 Effects of long range wakefields on a bunch 
train 

3.1 Energy deviations and kicks on the bunches 

The long range wakes due to HOMs can cause energy deviations and kicks on the 
bunches. A bunch train of bunches is shown in Figure [21 the notation for the 
offsets with respect to the reference axis of the accelerator is Xi and yi and the 
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direction of the longitudinal coordinate s (s = at the first bunch of the train) 
is also shown. It is assumed that all bunches have the same bunch charge q. In 
our investigation of the energy deviation and the kick on the bunch n within the 
bunch train we profited strongly from the analysis of long range wakes and beam 
loading by P. Wilson [6] . 



s 
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Figure 2: Schematic representation of a train of bunches with offsets Xi and 
with respect to the reference axis of the accelerator. 

The energy deviation of bunch n due to HOMs is a sum over all HOMs 
produced by all preceding bunches in the train. Additionally the the self-wake 
due to HOMs is also included {wI^\q) and wI^\q)), which is the HOM equivalent 
of the fundamental theorem of beam loading [6] . 

l^E{sn) = eqWl'^\0) (15) 

j<n 

+eq{xj + yn')wl'\0) 

+eqY^ {xn Xj + yn Vj) (s„ - sj) 

j<n 

+ ... 



As an example the energy deviation due the selfwake of a higher monopole 
mode with a frequency of / = 6 GHz is considered. The contribution to the 
energy deviation for a bunch with a population of = 6.5 x 10^ electrons is 
given by: 

1 /? 

= elnC 2 7r 6GHz - -. (16) 

2 Cj 
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An energy deviation of 196 eV corresponds to an impedance (R/Q) of 10 Q and 
a loss parameter of /cy^'' = 0.188 V/pC. The factor 1/2 is a consequence of the 
fundamental theorem of beam loading. 

Now consider the transverse long range wakefields. The kick on bunch n due 
to the dipole wake field is: 

On = ^^r^ i^j +yj ^yW±\^n- Sj) (17) 

^heam j<n 
+ ..., 

where En is the energy of bunch. The kick on bunch n due to one dipole mode 
with revolution frequency ui and damping constant ri is: 

On = 0n e^ + ^ ey)sin(a;i(s„-sj)/c) exp(-l/ri (s„-Sj)/c), (18) 

where the kick amplitude 6'„ on bunch n is defined as 

Qbunch ^ /inN 

^n = -^ c-7^ro, (19) 

with respect to an arbitrary reference offset tq. 



3.2 One dipole mode and a bunch train with constant 
offset 

It is instructive to consider the simplified situation of a bunch with a constant 
offset with respect to the "third harmonic" cavity. This corresponds to an injec- 
tion error or an misalignment of the cavity. Since the long range dipole wakefield 
is a linear superposition of HOMs it is sufficient to consider only one mode at a 
time in all analytical formulas. 

Modes from the first 3 dipole passbands with the highest values for R^^^ / Q are 
summarized in table [2], taken from reference [8], along with high R^"^^/Q modes 
of other azimuthal number. The kick amplitude 6 for the dipole modes has been 
calculated according to Equation assuming that all bunches have the same 
energy of 130 MeV, bunch charge of 1 nC and reference offset of 1 mm. 

Furthermore it is now assumed that the bunch to bunch distance At is con- 
stant: 

As 1 

At = — = n/fe— , (20) 

C Jfu 

where ffu = 3.9 GHz is the frequency of the fundamental mode and Ufh is the 
number of free buckets between bunches. The following bunch distances have to 
be considered for the operation of the injector linear accelerator: 
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//GHz 


m 


i?("V<5 / 


Gi / 

n 


V/(pCcm(2'^)) 


9 / yurad 


7.506 





23.3 


475.5 


0.55 




4.834 


1 


50.7 


277.3 


0.77 


1.22 


5.443 


1 


20.9 


426.2 


0.36 


0.50 


7.669 


1 


29.5 


470.7 


0.71 


0.71 


9.133 


2 


11.2 


402.9 


0.32 





Table 2: RF-parameters and kick amplitude of modes with high R^^^^Q. For 
the dipole modes (m = 1), the kick amplitude has been calculated for an beam 
energy of 130 MeV, a bunch charge of 1 nC and an reference offset (ro) of 1 mm. 



At / ns 


1 / At / MHz 


rifb 


200 


5.0 


780 


1000 


1.0 


3900 


2000 


0.5 


7800 


10000 


0.1 


39000 



Table 3: Typical bunch to bunch distance for the operation of injector linear 
accelerator 



The bunch distance can be translated into a phase distance S between bunches: 

S = uJiAt = 2TT^nfb, (21) 

Jfu 

where ui = 2ti fi is the frequency of the considered dipole mode. 

A small change in the dipole frequency due to fabrication tolerances will cause 
a large change in the bunch to bunch phase since the number of free buckets is 
relatively large. One obtains for a bunch distance of 1/At = IMHz: 

A5 = 27r-^n^bA/i (22) 

Jfu 

= 0.36°^ (23) 
kHz 

A change of 10° in the bunch to bunch phase after Ufb = 3900 free buckets 
corresponds to a frequency shift of about 20 kHz, which is much smaller than the 
expected variation due to manufacturing variations. 

A bunch to bunch damping constant for the kick voltage is defined as: 

"^ir^'^^'i""^- '''' 

where Qi is the Q-value of the dipole mode, which is usually dominated by the 
external Q-value of the HOM damper. We give most of our expressions both 
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in terms of d and of a = e"''. A plot of the damping constant as a function 
of the Q-value is shown in Figure [3] in a double logarithmic scale for the three 
modes considered in table [2] and a bunch to bunch spacing of 1/At = 1 MHz. 
With HOM dampers Q-values of about 10^ are achieved, which corresponds to 



Damping constant d for 1 MHz bunch spacing 





f = 4.8 GHz . 




f = 5.4 GHz 




f = 7.6 GHz 




d = 0.50 - 




d = 0.15 ■ 

















1000 10000 100000 le+06 le+07 le+08 le+09 

Q-value 

Figure 3: The damping constant d versus the Q-value for the three dipole modes 
from table [2l Additionally the damping constants d = 0.5 and d = 0.15 are 
indicated with parallel lines. 

a damping constant of about 0.15 (1/At = 1 MHz). If no HOM-dampers are 
mounted on the cavity the Q-value will be larger than 10^ and the damping 
constant d will be very small (~ 10~^). Some details of the damping constants 
for the modes from Table [2] are listed in the Table HI Using the above defined 





1/At = 1 MHz 




1/At = 5 MHz 


//GHz 




Q- Value 








Q- Value 






10^ 


5 X 10^ 


10^ 




10^ 


5 X 10^ 


10^ 


4.834 


1.5 


0.30 


0.15 




0.30 


0.06 


0.04 


5.443 


1.7 


0.34 


0.17 




0.34 


0.07 


0.03 


7.669 


2.4 


0.48 


0.24 




0.48 


0.10 


0.05 



Table 4: Damping constants d for three modes and different external Q-value and 
bunch-to-bunch distances of 1000 ns and 200 ns. 

parameters 6 and d, the energy deviation and the kick on bunch n G {1,2,3...} 
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due to one dipole mode are: 

AE^ = E(^^ + J2cos{5{n-j))exp{-d{n-j))j, (25) 

n-l 

On = e Y.MHn-j))exp{-d{n-j)), (26) 
i=i 

where 

E = -equJi-—ro , and 6* = — — —— Tq. (27) 
Q Eo/c Q 

To is a reference offset. The above expression fl25|) and fl26l) can be rewritten as 

AE„ = iQ + Re(^„; 

e„ = ^ Im(S'„), (28) 
with a sequence of complex sums Sn, defined as 

n—l n—l 

Sn = Y. exp ((n - j) /^) = E exp (j /^) . (29) 

i=i i=i 

The complex damping constant D is defined as D = i6 — d. The sequence S'yj 
may be calculated via a recurrence relation: 

Si = (30) 

or via an explicit expression for the sum of a geometric series: 

l-expiin-l)D) 1 

Sn = ^^TT > ^^TT for n^oo. 31) 

exp(— iJ) — 1 exp{—D) — 1 

Furthermore let i?„ be defined as: 

exp((n - 1)D) 
exp(— 19) — 1 

so that 

Sn = liin Sn — Rn- 

n— »oo 

The explicit expressions Re (/?„) and Im (i?„) are 

g-n.d (cos(n(5) — e"*^ cos((ri — 1) 5) ) 

Re(K) = ^ ^ — ^ (33) 

^ ^ 1 - 2e-'^cos(5) +e-2'^ ^ ^ 



e ""^ ( sin(n 5) — e sin((n — 1) 5) 

Im(i?n) = \ ' ' ' . (34) 

^ ' 1 - 2e-'^cos(5) +e-2'i ^ ^ 
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3.3 The functions Fji^n and F/„ 

The energy deviation AEn of bunch n and the kick 6n on bunch number n caused 
by the previous bunches are 

AEn = - Re (Sn)) = EFR,r,{6, d) 

and 

These equations define the functions Fji^n{S,d) and Fi^n{S,d), which depend on 
the bunch to bunch phase advance 6 and the damping constant d. A plot of the 
functions c?) and Fi^n{^,d) is shown in Figure H] for bunch number 10 and 

for a damping constant d = 0.15. 



Functions Fp „ and F, 




-180 -150 -120 -90 -60 -30 



30 60 90 120 150 180 



Figure 4: The functions d) and Fj^^i^, d) versus the bunch to bunch phase 

5 for d = 0.15 and n = 10 bunches. 



We have the following explicit expressions for these functions: 

— a" cos(n5) + a^~^^ cos((n — 1)S 



FR,n{S,d) 



Fi,ni5,d) 



(l-2acos(5) + a2) 

a sin(5) — a" sm{n6) + a"'~^^ sin((n — 1)5 
(1 -2acos(5) + a2) ' 



(35) 



(36) 



The bunch to bunch phase can be regarded as a random quantity since the 
interbunch phase difference S variation from expected variations in manufacturing 
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is much larger than ±27r. Consequently, it is valuable to calculate the average 
and the rms value of d) and Fi n{S, d) over the range — vr < 5 < tx. We are 

also interested in the absolute kick amplitude which is measured as the average 
of These are: 



{FR,n) 
{Fl,n) 
{\Fl,n\) 



1 

2^ 



d6FR,n{S,d) 



1 



1 7r J-TT 



(37) 
(38) 



TT JO 

1 1 - 



TT 



E 



A; 



exp(— A; d) 



k=\ 

1 / exp(— n (i) 



vr 



n 



i/(n, c?) +ln(coth(-)) 



(39) 



where the function i?(n, rf) is defined in terms of the hypergeometric function 
[9] as: 

H{n,d) = 2Fi{l,n;n + l;-exp{-d))- 2^1(1, n; n + 1; exp(-rf)) (40) 

Furthermore one obtains for the RMS-values 



1 exp(— nci) sinh((?7, — l)d) 
^ 4 ^ 2 sinh(d) 



1 + a2 - 2 a' 



2n 



\ 4(1 



(41) 



RMS(F,,„) = J^£jSFr,^{S,dy 



\ 



exp{—nd) sinh((r;, ~ l)d) 
2 sinh(d) 



(42) 



1 a2 - a2" 
^2(1^- 

For small (<^ ra) values of d we have: 



RMS(F,,„) 



In — 1 



for ^ 0. 



(43) 
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Because of the appearance of n — 1 in the square root here, it is instructive to 
plot RMS(F/ „,) as function of the bunch number minus 1. A plot of RMS(F7 „) 
for five different damping constants d in the range from 0.1 to 10~^ is show in 
Fig. [5] for bunch numbers n — 1 from 1 to 10^. The RMS- value of -F/^„ as a 



RMS(F,„(d)) 



100 





100 

bunch number n-1 



10000 



Figure 5: The functions RMS(F/ ,; 
d. 



versus n — 1 for different damping constants 



function of n quickly approaches the RMS-value of the asymptotic function Fi if 
the damping constant is relatively large. To be more precise we define the ratio 
r of the RMS- values as: 

RMS(Fj,„) 



Equation 



lim™RMS(F7,„) 
can be solved for the bunch n: 

ln(l — 



2 (n-1) ^ ^ 



n 



1 + 



2d 



(44) 



(45) 



where the bracket [ ] indicates the smallest integer not less than the expression in 
the bracket (ceiling function). In Tab. Owe have summarized the bunch numbers 
n for which RMS (-F/ ^) /RMS (F/) is larger or equal r = 0.95 using the damping 
constants from Fig. [5l For a given damping constant d we can find the bunch 
number n according to Eqn. fH5|) for which the RMS- value of asymptotic function 
Fj is a good approximation to the RMS- value of the function -Fj For a bunch 
train with 800 bunches the asymptotic function is only a useful approximation if 
the damping constant is larger or equal to 0.01 which corresponds to a Q-value 
which is smaller or equal than 10^. 
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d 


0.1 


0.01 


0.001 


0.0001 


0.00001 


n 


13 


118 


1165 


11641 


116397 



Table 5: Bunch number n for which the ratio of RMS values of the function F/^- 
and Fj is equal or larger than r = 0.95. 



3.4 The asymptotic functions and Fj 

For nearly all cases of interest, the functions Fji^ni^, d) and Fj^n{S, d) reach their 
large-n asymptotic limits Fji{6,d) and Fj{S,d) after only a small fraction of the 
bunch train has gone through the cavities. The expressions in this section are 
then useful. 

FR{6,d) = l+lim Re{Sn) 

1 ^-2d 



2(1 -2e-^cos(5) + e-2'i) 

sinh((i) 
2 (cosh(d) - cos((5)) 
l-a2 



2(1 -2acos(5) + a2 
Fi(6,d) = limlm(^„) 

71— >00 

e"'^ sin((5) 



(46) 



1 - 2 e-'^ cos((5) + e-2'i 

sin(5) 
2(cosh(rf) -cos(5))' 
a sm{6) 



— (47) 
(1 -2acos(5) + a2) ^ ' 

A plot of these functions is shown in Figure[6]for a damping constant d = 0.15. 

Fundamentally, we are exciting a simple harmonic oscillator with a train of 
(5-function pulses, and -Fr(5, d) is proportional to the response of the oscillator. It 
shows a characteristic resonance bell-shape, with the peak at the condition where 
the S pulses arrive at any sub-harmonic of the oscillator. The functions Fr((5, c?) 
and Fi{6,d) are asymptotic amplification factors along the bunch train; there is 
no bunch-to-bunch amplification of the energy deviation AE^ or kick 6n if these 
functions are smaller than one. 

The average and the rms values of -Fr((5, d) and Fi{S, d) for — vr < 5 < ti are: 

{Fr) = ^jy^FR{5,d) = \ (48) 
{Fi) = r d6Fj{6,d) = (49) 

2 TT J-TT 
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Functions Fj^ and Fj 




-4 I 1 1 ' 1 1 1 1 1 1 1 1 1 

-180 -150 -120 -90 -60 -30 30 60 90 120 150 180 

5 

Figure 6: The functions d) and Fj{S, d) versus the bunch to bunch phase 5 
for d = 0.15. 



(50) 



By rewriting the finite series for the average of the absolute kick function as 
a difference of two infinite series we obtain 



\Fi\) = - r d6Fii6,d) = -HcothA) 
n Jo 71 Z 



(51) 



1 2 
-ln(^) 

TT d 



Though the expectation (Fj) in the interval [— vr, tt] is zero, the expectation 
for the absolute value of the kick gets larger as d gets smaller. 
Furthermore one obtains for the asymptotic RMS-values 



RMS(Ffl) 



2tt 



dSFR{5, dy 



-^coth(rf) 



l + g^ 
^ 4(1 -a2) 



(52) 



^We have J^kLi ^exp(— fcrf) = -ln(l -cxp(— d)). 
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RMS(F/ 



27r 



d6Fj{6, dy 



-^(coth(rf) - 1) 



\2(1 



2v^' 

where the last approximation is vahd in the small d limit. 

The RMS of around its mean of 1/2 is equal to the RMS of Fj: 



(53) 



RMS(F«)2 - (1/2)^ 



RMS(Ff 



-^coth(d)-l 



(54) 



2Vd 



+ ■ ■ 



The RMS of the function Fr and Fi grows as l/\fd and will finally go to infinity 
if the damping is very small. 

Figure [7] shows the dependence of RMS{Fi) upon d. 



RMS(F,(d)) 



1000 




le-05 



0.0001 



0.001 



Figure 7: The RMS of the asymptotic Fi as a function of d. 
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In the no-damping limit, 



d—tO 



d 



2(1 -d) 1-6 



(55) 



which will go to d/6^ for small 6. Fj{6, d) approaches the function fi{6) which is 



Fi{6,d = 0) 

1 ,6, 
- cot(-) 

2 ^2' 



(56) 



This form is the result of taking the limit of Equation [3T] first for large n and 
then for small d. Taking just the limit for small d, 



o=^--M<n-m ^^^^^^ 



(57) 



exp(— «5) — 1 

is showing that at very small dampings the deflection, i.e. the imaginary part of 
Sn, will exhibit an oscillatory behavior before reaching equilibrium. 

Plots of function Fi{6,d) for different values of the damping constant d (0.1, 
0.15, 0.5) and the function fi{S) = Fj{6, d = 0) are shown in Figure [HI 
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Figure 8: The function Fj{S,d) versus the bunch to bunch phase S for different 
values of the damping constant d. 



The maximum of the function Fj is not aX S = like but at 

1 



(^max = arccosi 



~ h . . . 

cosh((i)^ 6 



(5^ 
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The maximum of the function at 6 = ^max is a function of the damping constant 
d: 

Max [Fj] (d) = FjiS,^,,, d) = / (59) 

2 smh (ct) 

Plots of functions RMS(F7)((i) and Max [F/] (rf) versus the damping constant d 
are shown in Fig. [91 
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Figure 9: The functions Max [F/] (rf) and RMS(F/-)((i) are plotted versus the 
damping constant d. 

The function -F/((5, d) is smaller than one for all possible phases 5 if the damp- 
ing constant d is equal or larger than 



= In Q (l + v^' 



0.481 . 



(60) 



This corresponds to a Q-value of about 3 x 10^ if a bunch to bunch spacing of 
1/ At = 1 MHz is considered. Such a low external Q-value is usually not achieved 
for all HOMs. But even if the damping constant is smaller than di the function 
F/((5, d) will be smaller than one for most phases (see Figure [H])- 
For d < di there exist two solutions of the equation 



Fi{5,d) = l. 

The larger one may be denoted as 5i. We have 

5i ~ arccos(3/5) — d^ 



(61) 



(62) 
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The probability P that Fi{5, d) is larger than one is therefore approximately 

Pfj>i ~ - arccos(3/5) ^ 0.3. (63) 

TT 

The amplification of the kick due to one HOM is expected to be smaller than one 
in about 70 % of the cavities of all possible random phases. Since several modes 
have to be considered the overall result will be different from this optimistic 
one-mode scenario (see the discussion in section [3^ . 

The probability that the function Fi{6, d) is larger than a constant A can be 
estimated from the proability that the function fi{6) is larger than A. In general 
we have: 

Pfi>a < Pfi>A = - arctan ( — j . (64) 

For many practical cases it is Pfi>a ~ Pfi>A if the constant A is smaller than 
Max[Ff]. 



3.5 The RMS kick over a bunch train due to one dipole 
mode 

In the previous section we have discussed the asymptotic kick on the bunch train. 
The average and RMS values of the asymptotic function Fr(5, d) and Fi{5, d) have 
been calculated for the phase 5. Now we drop the asymptotic limit, and average 
the kick On for finite bunch number n. We give the RMS kick over the bunch 
train as well. We retain the full form of the kick 



dnld = Fi{5, d)-\m{Rn) 



(65) 



where Fj and Rn are defined as in Equations HT] and [321 The average and the 
RMS of the kicks over N bunches in the train, normalized with respect to the 
kick amplitude 9 are 



1 ^ 

{6)16 = - - 5: Im (i?„) 

n=l 



(66) 



RMS{9)/9 



\ 



1 ^ / 1 ^ \ 

- E (Im {RnW - [j^ E iRn)j . (67) 



The analytic expression for the sums are: 



N 



T7 E (Rn] 
n=l 



-(Af+2) d 

2N 



1 + e-2'^ - 2e-'^ cos(5) 



cosh(c/) — cos((5) 



l)sin(5)- 



(68) 
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2 e'^ sin(iV(5) - e^'^ sin((iV + 1)5) - sin((A^ - 1 

1 N p-{N+2)d 1 

— y (Im(/?„))2 = ^. (69) 

iV^/ ^ N (1+6-2^-2 6-^008(5))' ^ ^ 

( - cosi5) + smHN d)] + 

\ sinn(aj J 



e 



-2{N+3)d 



2N l + e-2d_2e-'^ cos(5)2 
1 



l + e-id _ 2e-2rf cos(2(5) 

(2e(2^+i)'^(e2'^-l)cos(5) + 
g2ivd^g4d _ ^) ^^g(25) - cos(2 (iV - 1) 5) + 
e^"^ cos(2 {N + 1)6) + 2e'^ cos((2A^ - 1) 5) - 
2e^'^cos((2iV+ 1)5)' 



While the general analytic expression for the average and the rms kick over a 
bunch train of N bunches are rather complicated the expressions in the limit of 
no damping at all {d 0) are much simpler: 

Aveb{N,6) = \im{e)/e (70) 

1 fsm{N6) . \ 

— sm 



2 (cos(5) - 1) V N 



Rmsb{N,6) 



lim RMS{9)/9 

d— >0 



(71) 



2(1 -cos(5)) 



2sm{NSy 
Ar2 



+ 



sm(2 jV^) tan(<5)^ 
N 



;i -cos(5))^ 



Furthermore one may look for the limit of the functions Aveh{N, 5) and RmSblN, 5) 
for long bunch trains (A^ — > oo). Provided the limit is taken with constant 6, so 
that5A^ — > oo, 



lim Aveb{N,S) 

TV— >oo 



1 ,5, 
-cot(-) 

2 ^2' 

lim Fi{S, d) 

d — >0 



(72) 
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1 Jl-cos(5) 

lim Rmsb(N,6) = -^-^ ^. (73) 

2V2 (sin(f))' 

The limit of Aveb{N, 6) with respect to N is equal to the asymptotic amplification 
function Fj[6,d) in the limit of no damping, which diverges at 5 = 0. However, 
Aveb{N,S) is finite for all S > 1/N even if there is no damping. The functions 
Aveh{N, S) and lim^^o ^"1(6, d) are shown in Figure [TOl for a 100 bunch train. 




3.6 The RMS kick due to several dipole modes in an en- 
semble of cavities 

Based on the results of the previous subsections we want to calculate the expected 
kick due several modes in an ensemble of cavities. We want to find damping 
constants that limit the expected kick to less than the RMS beam divergence at 
the cavity. In this section, our estimates are based on the three dipole modes of 
Table [2j We assume that the electrical axis is the same for all modes and defines 
the axis of the cavity. The reference tq offset of the beam to the cavity axis is 
always 1 mm. Furthermore we assume that the betatron-phase advance between 
the cavities is small and the kicks from several cavities can be simply added. 



22 



3.6.1 An estimate based on the asymptotic RMS kick 

An estimate using an RMS approach and including a sum square combination of 
the different modes in the cavity is: 



Orms'= E r d5Fj{5,du 

, — , 1 ^ vr J -IT 



(74) 



fcG {modes} 

If the damping constants of aU modes are identical one obtains for 4 cavities: 



e 



RMS 



RMS(F,)(d) ■ ,AJ201 

k=l 

RMS{Fi){d) ■ 2.991 /irad 



(75) 
(76) 



The value 2.991/irad is based on the values in table [2l Equation fl76l) can be 
used as a criteria to determine the damping constant d. If we demand that the 
combined RMS kick should not exceed half of the single bunch divergence: 



2.991 /xrad ■ RMS(F/)(d) < lO/irad, 
one obtains for the required damping constant drmsw- 

drmslO = 0.0219 . 



(77) 



(7J 



The corresponding Q- values of the modes are listed in Table [6] and the function 
Fi{6,drmsio) is plotted in Figure [TT] for < 6 < ir. The RMS and maximum of 



//GHz 


Qext 


Qext 




{1/At = 1 MHz) 


{1/At = 5 MHz) 


4.834 


6.9 X 10^ 


1.4 X 10^ 


5.443 


7.7 X 10^ 


1.5 X 10^ 


7.669 


1.1 X 10^ 


2.1 X 10^ 



Table 6: Q- values of three modes for a damping constant d = 0.0219 and bunch- 
to-bunch distances of 1000 ns and 200 ns. 



the function Fj{6,d) for the damping constant drmsio = 0.0219 are: 
RMS{Fj){drmsio) = 3.33, Max [Fj] (d.^^io) = 22.9 . 
So in the worst but very unlikely case the total kick can be: 



tot 



Max [Fi] (drmsio) 
Max [Fi] (drmsio) 



k=l 

9.71 /irad = 222.4 /xrad. 



(79) 



^0) 
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The worst case assumes that the frequency of all 3 dipole modes in all 4 cavities 
are tuned in a way that the bunch-to-bunch phase is just equal to 5max where 
the function (irmsio) has it maximum. The probability that the function 
drmsw) is larger than the RMS value is: 



F7>3.33 



9.35' 



^1) 



Function Fj 




180 



Figure 11: The function -F/((5, d) versus the bunch to bunch phase 5 for d = 0.0219. 
The RMS value (3.33) of the function is also shown. The function Fi{6, d) is larger 
than the RMS value in the interval [0.09°, 16.9°] 



3.6.2 Alternative constraints on the asymptotic kick 

Instead of constraining the RMS kick one may put a constrain on the RMS value 
of the maximum possible kicks: 

Qrms^ = E ^ (Max [Fj] (4))'. (82) 

fc£{modcs} 

If the damping constants of all modes are identical one obtains for 4 cavities: 



e 



RMS 



Max [Fi] (d) 



\ 



k=l 



= Max [Fj] (d) ■ 2.991 /xrad 
< 10/irad 



(83) 



^4) 
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The RMS value of the maximum possible kicks is smaller than lO/xrad if d = 
dRMS = 0.15, since 

Max [Fi] (dRMs) < 3.34 . (85) 
The corresponding Q- values are summarized in Table [71 The RMS of the function 



//GHz 


Qext 

(1/At = 1 MHz) 


Qext 

{1/At = 5 MHz) 


4.834 
5.443 
7.669 


1.0 X 10^ 
1.14 X 10^ 
1.6 X 10^ 


2.0 X 10^ 
2.3 X 10^ 
3.2 X 10^ 



Table 7: Q- values of three modes for a damping constant d = 0.15 and bunch-to- 
bunch distances of 1000 ns and 200 ns. 

Fj{6,d) is 1.2 and therefore 

Orms = RMS{Fi){dRMs) ■ 2.991 fiiad 

= 3.59 fiiad , (86) 



9tot = Max [Fj] (dRMs) ■ 9.71 /xrad 

= 32.2 /irad. (87) 

Even in the worst case the total kick will only be a factor 1.5 larger than the 
single beam divergence of 20 /xrad. 

One of the most strict alternatives is to constrain the worst case of the total 
kick to 10 /irad: 

Otot < Max[F,](rf) ■ 4 ^^fc (88) 

k=l 

= Max [Fi] (d) -9.71 /irad 
< 10 /trad 

If the damping constant d is larger d^ax = 0.47 we obtain 

Max [F/] (rf) < 1.03 , (89) 

and the constrain Ofot < 10 /irad is met. The corresponding Q- values are listed in 
Tab. [HI For this very strict constraint we have: 

OjiMs = RMS{Fi){dmax) " 2.991 /irad = 1.69 /irad . (90) 
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//GHz 


Qext 

(1/At = 1 MHz) 


Qext 

{1/At = 5 MHz) 


4.834 
5.443 
7.669 


3.2 X 10^ 
3.6 X 10^ 
5.1 X 10^ 


6.5 X 10^ 
7.3 X 10^ 
1.0 X 10^ 



Table 8: Q-values of three modes for a damping constant d = 0.47 and bunch-to- 
bunch distances of 1000 ns and 200 ns. 

4 Monte Carlo Analysis 

The problem has also been attacked computationally with a simple Monte Carlo 
calculation. A worst-case analysis for a single cavity, and typical-case analyses for 
both FLASH and XFEL configurations have been made. The analysis provides, 
with certain assumptions, Qext requirements. 

The analysis does allow for the scatter of HOM frequencies that will no doubt 
result from variations within tolerances of the cavity dimensions and shapes that 
are part of the manufacturing tolerances. The analysis does not allow for the 
possibility that if HOM dampers fail utterly, fields may remain in the cavity for 
the relatively long time between bunch trains. Perhaps most importantly, there 
will be some active-feedback beam steering system that should ameliorate beam 
deflection effects and this remains still to be studied. 

4.1 Method 

The analysis is best explained by describing the objects from which it is made. A 
mode is basically a frequency, an value, an azimuthal quantum number and 
a decay time. A cavity is a collection of modes and the corresponding wakefield 
functions W\\ and W±. A beamline is a set of cavities, with the kicks and energy 
losses from each of the cavities summed. 

A technical note: the algorithm directly implements Equation [13] and for the 
bunch train lengths under consideration here, that requires taking the sin(6') and 
cos(^) for very large values of 6 indeed. Care has been taken to verify the numeric 
precision of these functions in this specific implementation of the algorithm. 

All the modes listed in Table [2] are included in the simulation, although the 
quadrupole mode has a negligible effect. The effect of the monopole mode on 
the beam energy is computed, although we have concentrated on the transverse 
dynamics. A sixth entry exists in the code for the 3.9 GHz operating mode, but 
it is switched off and is not used. 
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4.2 Single cavity worst-case analysis 

The worst-case is where a cavity has a HOM with a peak in Fi that is directly 
on beam resonance and the external dampers have failed so that damping is 
provided by the intrinsic power dissipation of the Nb surface alone. In this case, 
the asymptotic limit is not reached; the bunch train is not long relative to the 
Wakefield time scale. The surface resistance is modeled from 3.9 GHz with an 
Rbcs of 20nf2 scaled with the square of the frequency ratio, plus the residual 
resistance, i?0; of 2{]nVL. The Q-value is then obtained by dividing Gi of the 
mode by this resistance. 

The monopole mode will, with over 800 bunches in the FLASH configura- 
tion, take about 800 kel^ out of the last bunch in the train. The maximum of 
Fj corresponds to about 468 Yiz from the exact on-resonance condition for all of 
the dipole modes 0. The deflection from the dipole mode of highest -^^^ in this 
case is shown in Figure [T21 the other modes have the same shape but with the 
vertical scale proportional to Figure fT3l shows the crabbing angle, defined 

as the deflection angle (in the lab frame) at the l(Jz bunch head minus the de- 
flecting angle at the bunch center. This is computed by changing, in effect, the 
trailing distance of the witness by one bunch-length, and taking the difference in 
displacement from the centroid displacement. 

We do not have a clear understanding as to how this can effect the lasing 
process, but note that the angles involved are typically smaller than the deflecting 
angles. The energy loss in this dipole mode at this pseudo-resonance condition 
is only bkeV . 

It is clear that even a single cavity hitting the worst-case will cause lasing to 
stop. 

4.3 Typical-case analysis, FLASH beamline 

In the typical-case analysis we use the Monte Carlo method to examine the prob- 
ability distribution for beamline performance. 

A virtual beamline is constructed of 4 cavities. The frequencies of the modes 
are shifted from the nominal simulation results of Table [2]by drawing on a uniform 
random distribution. The width of the distribution spans the full range of 5, 
corresponding to a frequency shift of —0.5 MHz to -1-0.5 MHz in the FLASH 
case, and l/5th that for the XFEL injector. While we do not have enough 
cavities to really check, this is thought to be about the scale on which the scatter 
will actually be. 

Simulated results such as in Figures [12] and [131 along with the corresponding 
energy loss plot are recorded. Then a new beamline, constructed with new calls to 

^In the long bunch train hmit, the maximum would occur at A/ = f /2Q. However, without 
damping, 800 bunches is not near the long bunch train limit. Additionally, there are large-N 
oscillation effects, as described in section 
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Figure 12: Bunch deflection as a function of bunch number, worst case analysis 
of the 4.834 GHz mode. For this mode, Q = 5A x 10^ d = 2.77 x 10-^ and 
^=1.22 mrad. 
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Figure 13: Bunch crabbing as a function of bunch number, worst case analysis of 
the 4.834 GHz mode. 
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the pseudo-random number generator is created and again, variation of deflection, 
crabbing, and energy loss with respect to bunch number is determined. The 
process is repeated 5000 or in some cases, 10000 times. At each bunch number, 
we have an average and an RMS, over those beamlines, of deflection, crabbing, 
and energy loss. 

The average deflection over an ensemble is easily seen to be zero, as deflection 
to the left is just as likely as deflection to the right. The average energy loss is also 
zero; the sin of Equation [T3b is replaced by the cos of Equation [T3k. The widths 
of the distributions of deflection angle, crabbing angle and energy loss describe 
at the la level what kind of performance one can expect when the machine is 
actually turned on. 

Clearly, a design that will keep deflection down to our ±10/irad goal at the 
68.27% (Icr) confidence level is risky; we need the RMS to be well below the goal. 
How far below is a difficult decision involving overall program risk. Notwith- 
standing, a statistical analysis is able to provide us with a good sense of what 
levels of external damping we need to have. 

Again, we do not obtain an acceptable result on just cavity self-damping alone. 
If the HOM design fails across the board, deflections of ~ ±50yurad appear at 
the la level by the end of an 800 bunch train in the 4 cavity, 5 mode FLASH 
model, as shown in Figure [TH Next, we determine how far we smust lower Qext 
to stabilize the beam. 




100 200 300 400 500 fiOO 700 SOO 

Bunch 



Figure 14: The one-sigma deflection profile as a function of bunch number, aver- 
aged over ensemble of 5000 virtual FLASH beamlines without external dampers. 

For the 7.506 MHz monopole mode, lowering Qext to a bit below 1 x 10^ causes 
the RMS of AE to flatten out at large bunch numbers to about 23 keV. 



29 



For the dipole modes, there is some allocation of the ilO/xrad target into 
the three modes; we have to also allocate fractions of the deflection budget to 
the different deflecting modes. If we require that the RMS deflection flattens 
out to a level where zblO/irad corresponds to 3a of the total deflection, and then 
allocate the deflection budget evenly among our three large dipole modes, then 
each mode must contribute 10/irad/(3y(3)) = 1.9/irad in the asymptotic limit 
for the 4 cavities, as the contributions from the included modes are added in 
quadrature. 

The Monte Carlo calculation differs from the analytic form of section 13.6.11 in 
that 

1. In the MC method, the asymptotic limit is not assumed. It will in most 
cases be justified because the damping is strong. 

2. In the MC method, angles are summed over modes and cavities, and then 
an RMS is taken; in the analytic form, the sequence of these two operations 
is inverted. 

3. In the analytic form, the damping constants of all the modes are forced to 
be equal; in the MC method, the deflection in the long-bunchtrain limit are 
taken to be equal. 

4. Equation [77] sets lO/irad to 1 times the RMS deflection; section 13.6.21 dis- 
cusses alternate, tighter constraints. The MC method uses a 3a constraint. 

The Monte Carlo results can be reproduced with analytically. With the re- 
quired values of Qext, the asymptotic form of Equation [531 permits direct solution 
for Qext using the definition of 6 and d. These relations lead to Qext for the 3 
modes of 3.7 x 10^, 2.5 x 10^ and 1.7 x 10^ or about 70-80% of the MC values. 

The Qext values that produce deflections of the specified level in the Monte 
Carlo, found basically by trial and error, are shown in Table [9], along with the 
results of section [3.6.11 when a 3a constraint is used, the results of Phillipe Plot's 
(unpublished) calculation and the recent measurements from prototypes (TU]. The 
system is well damped in the simulation, as shown in Figure [151 



4.4 Typical-case analysis, XFEL injector beamline 

The typical-case analysis for our canonical XFEL configuration shows, as ex- 
pected, that the RMS deflection due to HOMs varies inversely as the beam energy 
and approximately as the square root of the number of cavities. This makes the 
deflection about 2/3 of what it would be in FLASH at the same bunch spacing. 

Moving the bunch frequency up to 5 MHz changes the relative contribution 
of the dipole modes by changing d and thereby increases the accumulated deflec- 
tions. Where the asymptotic RMS values of the deflections, scaled to 24 cavities 
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Frequency 


Azimuthal 




Qext 






(GHz) 


number 


MC method 


3a Analytic 


Piot spec. 


Measured 


7.506 





1.0 X 10^ 






1 X 10^ 


4.834 


1 


5.4 X 10^ 


9.0 X 10^ 


2.8 X 10^ 


1 X 10^ 


5.443 


1 


2.9 X 10^ 


1.0 X 10^ 


1.0 X 10^ 


5 X 10^ 


7.668 


1 


2.1 X 10^ 


1.4 X 10^ 


2.5 X 10^ 


1 X 10^ 



Table 9: Qext requirements from this Monte Carlo study, the analytic forms of 
13.6.11 when calculated for a 3a constraint, Phillipe Plot's analysis, and measure- 
ment from prototypes. The FLASH beam parameters are used. The value of 
drmsio/3 is 0.169. 



■D 

B.003 
0.002 
0.001 


-0.001 
-0.002 
-0.003 

100 200 300 400 500 600 700 800 

Bunch 



Figure 15: The one sigma deflection profile as a function of bunch number, av- 
eraged over ensemble of 10000 virtual FLASH beamlines, each consisting of 4 
cavities and 5 well-damped modes. 

in a 500 MeV beam would otherwise be 1.2/xrad, the introduction of 5 MHz bunch 
frequency results in deflections of about 2.8/irad. To restore 3a lOyurad perfor- 
mance, we need the Qext values of table [101 

The overall deflection profile still shows a well-damped system; the asymptotic 
one sigma deflection of 3.39/irad is reached within 60 bunches. The asymptotic 
energy change is ~ 58keV. 
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Frequency 


Azimuthal 


Q 




(GHz) 


number 


MC method 


3a Analytic 


7.506 





1.0 X 10^ 




4.834 


1 


2.4 X 10^ 


1.8 X 10^ 


5.443 


1 


1.3 X 10^ 


2.0 X 10^ 


7.668 


1 


1.0 X 10^ 


2.8 X 10^ 



Table 10: Qext requirements from this Monte Carlo study and the analytic forms 
of 13.6.11 when calculated for a 3a constraint. The XFEL beam parameters are 
used. 

5 Conclusion 

We have studied the Qext requirements for the "third harmonic" cavities using sets 
of beam parameters typical of FLASH and XFEL operation. A key assumption 
is that lasing will cease when the deflection due to wakefields approaches the 
divergence of the beam at the cavity location. We have taken the divergence 
to be ±10//rad for both beam parameter sets. The XFEL optics in the vicinity 
of the "third harmonic" cavities is not finalized as we write, but the ilOyurad 
condition is of the correct scale. 

The reader need also be aware of the 'program risk' issue: what level of 
statistical confidence that this ±10/irad specification be met without component 
replacement or change of beam parameters is required? Here, we have selected 
a 3a level of confidence. One might choose to argue that the 3a choice is too 
conservative. Choosing a la requirement relaxes the damping requirements by 
about an order of magnitude. 

The results are based basically on three modes of high beam-cavity coupling. 
Adding in quadrature a number of other modes with lower R^^^Q has little 
influence (~25%) on the damping required. 

Both analytic and Monte-Carlo based analysis have been done. Both allow 
for manufacturing defects, but do not allow for the action of any kind of active 
beam steering system and are quite conservative in that regard. 

The results of the two analyses are broadly consistent and are summarized in 
tables M and [TUJ It is encouraging that the Qext values measured on prototypes 
are better than the required values. 

For damping values of Qext in the 10^ range, which are typical for functioning 
HOM mode dampers, the deflections reach their asymptotic values after a few 
lO's of bunches. In this case, our analytic results are particularly easy to use, 
and we summarize them here. 

The angular kick on bunch n G {1, 2, 3 . . .} due to a single mode of angular 
frequency a; and quality number Q in a train with bunches At seconds apart is 
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On = eFi{S,d) where 



7^ 6 Qbunch ^ 

On = — ^-?r''o 

6 = iuAt 

The function F/, which describes the bunch-to-bunch amphfication of the 
wakes, has a maximum at 



dmax = arccos( — — — -) ^ d- — + . . . 

cosh(aj 



2 sinh(c/) 
and an RMS of 



RMS(F,) = ^l-£j6Fi{6,d) 

= y{coth{d) - 1). 
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